Vacuum entanglement enhancement by a weak gravitational field 



M. Cliche and A. Kempf 

Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1 

(Dated: August 31, 2010) 

Separate regions in space are generally entangled, even in the vacuum state. It is known that this 
entanglement can be swapped to separated Unruh-DeWitt detectors, i.e., that the vacuum can serve 
as a source of entanglement. Here, we demonstrate that, in the presence of curvature, the amount 
of entanglement that Unruh-DeWitt detectors can extract from the vacuum can be increased. 

PACS numbers: 03.67.Bg, 03.70.+k 
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I. INTRODUCTION 

Two Unruh-DeWitt detectors that interact with a 
quantum field in the vacuum state have access to a renew- 
able source of entanglement [IH5], namely by swapping 
entanglement from the quantum field. In this context, it 
was recently shown [B] that, in an expanding space-time, 
the entanglement of the vacuum decreases significantly 
due to the effects of the Gibbons-Hawking temperature 
[7j. While this example showed that gravity is able to 
act as a decohering agent, we will here show that gravity 
can also act to enhance entangelement-related phenom- 
ena. Namely, we will show that a weak gravitational 
field, such as that caused by a planet, can enhance the 
extraction of entanglement from the vacuum. 

This article is organized as follows. In Sec. |H]we re- 
view the extraction of entanglement with Unruh-DeWitt 
detectors in Minkowski space-time 



III 



we review 

Wl 



In Sec 

the Newtonian limit of general relativity and in Sec 
we look at Unruh-DeWitt detectors in the presence of 
weak gravity. Then, in Sec. |V]we compute the first or- 
der correction to the propagator on the perturbed back- 
ground. In Sec. |VI| we calculate explicitly the entangle- 
ment between the two detectors near a spherically sym- 
metric star. In the last section, we propose extensions. 

We work with the natural units H = c = G = 1 and 
the Minkowski metric rj^ = diag(—l, 1,1,1). We denote 
the coordinate time by x° — t while the proper time is 
denoted by r. Wherever necessary to avoid ambiguity 
we will denote operators O or states which live in the 
Hilbert space of the j'th subsystem by a superscript 
(j), for example, and \tp^'). Orders in perturbation 
theory will be denoted by a subscript (j), as in, e.g., P = 
P(o) + P(i) + 0(e 2 ). We work in the interaction picture. 



II. VACUUM ENTANGLEMENT 

Let us first briefly review vacuum entanglement with 
Unruh-DeWitt detectors [TJ |21 E] ■ To begin, let us denote 
the overall Hilbert space by U = % (1) <g>'H (2) <g)'H (3) , where 
the first two Hilbert spaces belong to two Unruh-DeWitt 
detectors and where the third Hilbert space is that of a 
quantum massless scalar field. The total Hamiltonian H 



of the system with respect to the coordinate time t is 
H = Hp + Hp, + Hi n t 



H F = 



1 



d 6 x 



n 2 (x) + (V(£(a;)r 



Hr 



^2[(E g + AE)\e^)(e^\ 



dt 

2 

H int (t) = ^Mrj(t))(\e (j) )(9 (i) \e iAE ^ 

3=1 

+ | fl W)><eW|e-* A *-iW)^(a !j (r J (t)))*# 



dt 
(1) 

where Hp is the Hamiltonian of a free massless scalar 
field, Hp> is the Hamiltonian of the two detectors, 
Hint(t) is the interaction Hamiltonian [5] in the in- 
teraction picture, otj is the coupling constant of the 
j'th detector (j G {1,2}), 4>(%j(Tj)) is the field 
at the point of the jth detector and m^'(Tj) :— 
(\e^){g^\e iAET i + \g^)(e ( ^\e- lAE ^) is the monopole 
matrix of the jth detector. The function T)(jj) will be 
used to describe the continuous switching on and off of 
the detectors and Tj is the proper time of the jth detec- 
tor. 

Let us first consider the special case where T\{t) = 
T 2(t) — r(t) such that the evolution operator U = 
rp e —ifdtH int (t) ac ting on states takes the form: 



U = Texp^ -i dr 



air?(r)m (1) (T)^(xi(r)) 



+a 2 r] (t) (t) (j) (x 2 (r)) 



(2) 



We assume that the initial state of the system is 
|0g( 1 )g( 2 )). After the unitary evolution of the total sys- 
tem, we trace out the field and obtain at 0(a 2 ) [TJ [B] : 



P f 2) = Tr {3) (U\0gMgW)(09 {1) 9 {2) \tf 



X 

Pi Y 

Y' P 2 

X* I-P1-P2 



+ Q(a 4 ) 



(3) 
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in the basis |e«e( 2 )), \eWgW), |g (1) e( 2 )) and | 5 (1) <? (2) ) 
The matrix elements Pi , X and Y read 



Pj = a. 



31 

oo />oo 



2 ; dr I dT'v(T)v(T')e- iAE ^- T '^ 

J — oo J — oo 

xD{x 3 {t),x 3 {t')) (4) 



X = -aian I dr dr' r]{T)T](T')e iAE ( T+T ") 

J —oo J —oo 

x(D(x 1 (t),x 2 (t'))+D(x 2 (t),x 1 (t'))) (5) 

/OO t* oo 

dr / dr'77(T)7 ? (T')e- lAjE(r - T ' ) 
-OO J — OO 

xB(h(t),i 2 (t')) (6) 
where D(x,y) = (0\cf)(x)<fi(y)\0) . To measure the entan- 

f 1 2) 

glement of ' , we use the negativity [§] which gives: 

N = max (J(Pl - P 2 ) 2 + 4\X\ 2 -P 1 - P 2 , o) 

+0(« 4 ). (7) 

In order to obtain more explicit results, let us consider, 
for example, the case where the detectors are inertial and 
separated by a constant proper distance L p in Minkowski 
space-time, in which case x°(r) = r. Let us also assume 
that «i = a 2 — ol such that in Minkowski space we 
have Pi = P 2 = P and N = 2max(|X| - P) + 0(a 4 ). 
For simplicity we choose the switching functions to be 
gaussian: t)(t) = e~ T /( 2<T "). 

We will call X the exchange term and Pj the local noise 
term. This is because X can be interpreted as describing 
the exchange of virtual quanta between the two detectors, 
and Pj can be interpreted as describing the detection 
of virtual quanta by detector j. In order to allow the 
introduction of gravity (which will enter mostly through 
the propagator), it will be useful to view Pj and X as 
functions of the propagator. To this end, let us already 
ensure that the time ordering is respected. For X this is 
straightforward since the time integrations respect time 
ordering by construction. We can still simplify X by 
using the variable change s — r — t' and u = t + t' such 
that we have 



X 



/ dse-* 2 /^ 2 



x (G(fi,£ 2 ,s) + G(x 2 ,xi,s)) 



(8) 



where G(x,y) = (O|T0(x)^(y)|O) = G(x,y,x° - y°) is 
Feynman propagator [10] . For Pj , we introduce a conve- 
nient change of variables for the double integral over the 
(t, t') plane |llj . making u = r, s = t — t' in the lower 
half-plane r' < r and u = r', s = r 1 — r in the upper 
half-plane r < r', Pj becomes: 

Pj = 2a 2 a^ (^J™ dse- s 'K i ' T *>- i * Ea G(x j ,x j ,s?j . 

(9) 



In Minkowski space-time we use the Boulware vacuum 
such that the propagator is given by 



G(x,y) 



-1 



4tt 2 [(a; - y°) 



\x-y[ 



(10) 



where lim e _j.g+ is implicit. Using Eq. (10 1 in Eq. ([9| and 
Q we obtain the local noise and the exchange term in 
Minkowski space-time [B], 



47T 



X = 



AE^aerfc(AE<7) 



4L r , 



r A B v-i;/4^ erfc 



2(7 



(11) 

(12) 



where L p = \x\ — x 2 \ and erfc(a;) = 1 — erf(x). In 
the regime AEa — > oo and L p /a — > oo, we have P rs 

"sZaew and \ X \ ~ ^^Lli " ■ Therefore, to get a 
non vanishing negativity, we need AEa > To opti- 
mize the negativity in that regime, we set d ^ - = 0, 

which yields AE opt w ^ (l + 2<j 2 /L 2 p ). The resulting 
negativity is then N "~ '" ,T 



* opt 



,.2„4 e -i^/4o 



III. NEWTONIAN LIMIT 

Let us now briefly review the Newtonian limit of gen- 
eral relativity, see e.g. [T2]. In this limit we can write the 
metric as g a p = r) a p + h a p where \h a p\ <C 1. Note that 



under a small change of coordinates x p 



^ the 



term h a p has a gauge transformation h a p — » h a p — dp^ a — 
d^p. Let us define the quantity ^ := /i^ - Tj^h^/2. 
To simplify the Einstein equation, we choose to work in 
the Lorentz gauge in which hf tb = 0. In this gauge, the 
linearized Einstein equation reads d a d a h^ v — — 167rT M1/ . 
In the Newtonian limit the gravitational field is too 
weak to produce velocities near the speed of light, thus 
only the T 00 component of the stress-energy tensor con- 
tributes significantly and we can make the approximation 
d a d a fa V 2 . This means that the Einstein equation can 
be approximated as d a d a h m sa V 2 ft°° « —Wirp. From 
this we conclude that the dominant component of is 
h 00 , such that in terms of h a/B we have h 00 = h u = h ot) /2. 
Thus, the line element takes the form: 



ds 2 



-(l-h 00 /2)dt 2 

+ (l + h 00 /2) (dx 2 + dy 2 + dz 2 ) . (13) 



Now assume we have a compact object, say a star of 
dark matter that does not interact with the quantum 
field and is of radius R and of constant density p = 
3M/(4irRl). We solve V 2 /) 00 w -16np with the usual 

boundary conditions /t 00 (|r| — > oo) = 0, dh g ^ - (r = 0) = 
and with the continuity conditions /i 00 (|F| — > R a — e) = 
h 00 (\r\ -> R +e), ^(\r\ fl -e) = i? 0+e ) 



3 



(14) 



in the limit £4 0. This gives 

£00 fi * = when|r1<i? , 
' \ ff when\r\>R 

so to have \h a p\ <C 1 we require M/R <C 1. 



IV. DETECTORS ON THE CURVED 
BACKGROUND 

Let us now consider the two Unruh-DeWitt detectors 
on the background of the weak gravitational field. We 
assume that the two detectors and the center of the star 
are all on a same axis. Therefore, detector 1 is located 
at a fixed distance r\ from the center of the star and 
similarly detector 2 is located at r 2 = r\ + L from the 
center of the star. This means that their proper times do 
not coincide 7i(t) ^ t 2 (£), so we may write the evolution 
operator as 



U 



Texp< 



d,T\a 



i7(n)mW(ri)0(a;i(ri)) 



-V (r 2 (n)) mP> (r 2 ( n )) 4> (x 2 (r 2 ( n ))) 



(15) 



and using Eq. ( 13 ) we have 



'l-ji°°(ra)/2 



= n 1 - 



l-/i 00 (r 1 )/2 

/i 00 (r 2 ) E 00 (n) 



O([/i 



(16) 



To simplify our analysis we want to avoid this blueshift 
effect. To do this, we assume that the two detectors are 
close enough such that their internal clocks have the same 
speed at first order in perturbation theory. This will be 
so if \h 00 (r 2 )/4- h 00 (n)/4\ < O([h 00 (r 2 )} 2 ) which for 
detectors outside the star gives L < 16M. Under that 
assumption, we have r 2 = Ti (1 + ([h 00 ] 2 )) such that 
one can easily verify that Eq. Q and Eq. ftty still hold 

uptoOp] 2 ). 

We are therefore left with two first order contributions 
to the exchange term X and the local noise term Pj , the 
first one which we denote by Xm\ and Pj(i) is essentially 
a result of the time dilation caused by the star and the 
second one which we denote by SX^ and SPjm comes 
from a modification of the propagator on the curved back- 
ground. Let us denote the perturbative expansion of the 
propagator as G(x, y) = G^(x, y) + G(i)(x, y) and since 
it is widely believed that a Boulware-like vacuum is the 
right vacuum for a quantum field in a Newtonian gravi- 
tational potential [121 [H], we use Eq. (10) for G( )(a:,y). 



We can easily evaluate the contributions Xn\ and Pjm 
by first noting that 



, l-h 00 (r j )/2 



(17) 



Thus, when L < 16M we have using Eq. ( 10 ) 



G m (\x 1 (r)-x 2 (r%x° 1 (r)-x 2 (r')) = 

(1 - ft 00 (n)/2) G (0) (L p (1 - h ao ( n )/2) ,r-r') 

+O([h 00 ] 2 ) (18) 



where 



Lp : — 



l + h 00 (r)/2dr 



L(l + h 00 ( ri )/A) 



(19) 



is the proper distance between the two detectors. Hence, 
when we put this back in Eq. Q and Eq. Q we have 
the first order corrections 



F m = o — ^(o) 



(20) 



A(i) = — IX (0 )+Lp^— ] \ 1 ) 



where the zeroth order terms are given by Eq. (11 1 and 
p). 



V. CORRECTION TO THE PROPAGATOR 

In this section we compute the first order correction 
to the propagator on the perturbed background. The 
first steps of our calculation can be found in |13j . To 
focus on the correction caused by gravity, we assume that 
the field is minimally coupled to curvature and to the 
matter that composes the star. Under these assumptions, 
the propagator is a Green's function of the Klein-Gordon 
operator 



O x G(x,y) 



i5(x, y) 



(22) 



where D x f(x) = -^d^[^gg^ v d v f{x)]. The first or- 
der correction to g is g = — 1 — = —1 — h 00 . Using 
G(x,y) = G( ) (x,y)+ G (1) (x, y) we have: 



=d u 



y/T+h%(rf u -hi iV ) 



x3„ (G (0) (x,y) + G {1) (x,y)) 



V 1 + K 



(23) 



Expanding everything to first order only and using the 
fact that G(o) (x, y) solves the zeroth order equation, we 
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obtain 

- h^d^G^x^) + n {Q)x G {1) (x,y) 

-d^d„G {0) (x,y) + d^/2)r]^d u G i0) (x,y) 

= -i5(x,y)h«/2 (24) 

where we used 0(0)1; := -q^d^d^. Using again the fact 
that i5(x,y) = U\(Q) x G(p)(x,y) we can simplify the previ- 
ous equation, 



O( )xG {1) (x,y) 



= d^ v d u G (0) (x,y) 

+h^d tl d v G {0) {x,y) 
= h^d„d u G {0) (x,y) 
= h 00 (x)d 2 x0 G {0) (x,y) 



(25) 



where we used the fact that we are in the Lorentz gauge 
and that in the Newtonian limit h 00 is the dominant 
component of . Note that since the space-time we 
consider is static and asymptotically flat, the propagator 
G(x, y) can be seen as the analytic continuation of the 
unique Green's function on the positive definite section 
[13] . Since this holds order by order in perturbation the- 
ory, at first order perturbation we can use G( ) as the 



inverse of D(o) such that: 



G (1) (x,y) = -i / d*zG {0) (x,z)h w (z)di G (0) (z,y).(26) 



This equation gives us explicitly the first order correction 
to the propagator. It is clear from this equation that the 
entire space-time perturbation will modify the propaga- 
tor, and the most significant contribution will come from 
the patch of spac e-tim e near x and y. We now insert 
G(o)(x,y) in Eq. (26 1 and using the fact that h 00 (x) is 



independent of time, we obtain 



G(i){x,y) 



16tt 4 



dz°d 3 zh 00 (z) 



H(z + s)' 

(z - zi) 3 (S - z 2 ) 3 (z ~ z )(z + z ) 



(z - zi) 2 (z - z 2 ) 2 (z ~ z )(z + z ) 



(27) 



where we use the definitions Z x := \x — z\, Z y := \y — z\, 
s := x a — y , z := z° — x , z a := X + ie, z\ := —s + Y + ie 
and z 2 '■= — s — Y — ie. We can then perform the z° in- 
tegration with the residue theorem. We choose a closed 
contour in the upper half of the complex plane and the 
upper part of the contour is equal to zero because the in- 
tegrand vanishes sufficiently rapidly as z° = Re l ° 

R— foo 



We thus have: 
G(i)(x,y) = 



1 

8^ 



d 3 zh 00 (z) 



8(z + s) 2 
(z - Zi) 3 (z - z 2 ) 3 2z 

(S + s) 2 



4-4— 
dz 2 



(z - z 2 ) 2 (z - z )(z + z ) 
2 



(z - zi) 2 {z - z 2 ) 2 2z 

d ( 2 

dz \(z - z 2 ) 2 (z - z ){z + z a ) 



1 

8^ 



d d zh m (\z\) 



3(s 2 + Z x Z y )(Z x + Z v ) + Z 3 + Z 3 
{Z X Z V + ie){s 2 - [Z x + Z y + ie] 2 ) 3 

(28) 

Using the above equation we may now evaluate 5Pj(i) 
and 5X^. For SPj^, we have Z x = Z y = Z, such that 
the correction to the propagator can be greatly simplified 
with a simple change of variables 



G ( i)(f,f,s) = 



1 



2r7r 2 



3s 2 



i rr+R 

dRRh 00 (R) / dv 

J\r~R\ 

f Av 2 



(s 2 — Av 2 — it) 3 



(29) 



where r is the distance between x and the center of 
the star and s = x° — y°. The v integral can be per- 
formed analytically. Note that Eq. Q and Q were 
derived for detectors in Minkowski spacetime, where 
x°(r) = y°(r) = r, but since we are only interested at the 
first order perturbation, the effect of the time dilation in 
the corrected propagator would be a second order term 
which we neglect. We can thus use Eq. ^ and ^ with 
the first order correction of the propagator and with no 
time dilation, that is x°(t) — y°(r) = t. For the same 
reason, we can also use at this order L p = L. Therefore, 
we can put Eq. ( 29 1 in Eq. ([9| and we obtain the first 
order correction to the local noise 6Pjny. 




dse- s2 /^-* AEs 



2{rj +R) + (s - ie) 
2{ rj + R)-(s- ie) 

R\ + (s-k)\\ 1 



dR 



-R\-( 
2( rj +R)(2{r 




(s - ie) 2 (4(fj + R) 2 - (s - ie) 1 

2\r ] -R\{2{r. J -R) 2 -s 2 ) 
(s - ie) 2 (4(rj - R) 2 - (s - ief 



..(30) 
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Similarly for the exchange term SX^, we put Eq. (28) 
in Eq. ^ and we then use a simple change of variables 
to obtain 



SX (1 



a 2 aJ^e-° AE r°° 



(i) 



27r 2 ri 



II 



oo rr\-\-R 

< / dRRh 00 (R) / dvivi 
la J\n-R\ 



X 



3(s 2 + viv 2 ){vi + v 2 ) + v{ + vl 
(viv 2 + ie)(s 2 — [«x + v 2 + ie] 2 ) 3 



where v 2 



(31) 



The 



s integration can be performed analytically, such that 
we are left with a relatively simple expression for SX^ 
which involves only two integrations: 



SX {1) = 



a <7\/7re ' 



27r 2 ri 



n+i? 

dRRh m (R) I dvi 
\n-R\ 



rerfc 



VI. 



I6(v 2 + ie)a 4 1 V 2cr 

xe -K+^)) 2 /(4. 2 ) [2a 2 -( Ul+W2 ) 2 ] 

-2V^a(v!+i; 2 )l. (32) 



NEGATIVITY ON THE PERTURBED 
BACKGROUND 



We now have all the tools to compute explicitly the cor- 
rected negativity. Using the /i 00 (|r|) of Eq. (14) in Eq. 



(|30|) and (|32j), we can find <5-P,(i) and SXn) by numeri- 



cally evaluating the remaining integrals. Pjm and Xn 
can then be evaluated exactly using Eq. (20) and (21 
such we end up with the full noise term Pj = P(q) +Pjn\ H 
SPj(i) and the full exchange term X = X(q) +Xm+6Xri 



using Eq. (Ill and 112 



(i) 

for the zeroth order terms. This 
le negativity between the two de- 



allows us to compute t 
tectors using Eq. 

Numerical evaluations indicate that \X\ linearly in- 
creases with the strength of the gravitational potential 
M /R a of the star while Pj linearly decreases with M/R . 
Therefore, the negativity N linearly increases with the 
strength of the gravitational field M/R a . In a similar 
fashion, numerical evaluation of |X| and Pj indicate that 
the correction to the negativity N decreases roughly like 
Ro/r\ as ?'i/i? — > oo but remains positive, see Fig. (JlJ. 
On Fig. we chose two sets of parameters L p /a and 
aAE. With the first one, the negativity N is non-zero 
even without the gravitational field such that we only 
have entanglement enhancement by gravity. With the 
second set of parameters, we have N = and |X| « P 
without the gravitational field such that in that particu- 
lar case we not only have entanglement enhancement by 
gravity but also entanglement creation by gravity. 



1.5 



N' 1 



0.5 



10 



15 



n/Ro 



FIG. 1: N' = 8ir 2 N/a 2 as a function of r 1 /R . Here aAE = 
0.00674, AER = 1 and M/R a = 0.001. The upper (red) 
curve correspond to L p /R — 0.0095 and the lower (blue) 
curve to L p /R = 0.01. The upper and lower dashed lines are 
the asymptotes r 1 /R a — > oo. 



We may heuristically interpret this phenomenon by 
looking at the local noise term Pj and the exchange term 
\X\ separately. Since the gravitational field increases the 
momentum of virtual particles near the star (as seen by 
the fixed detector j), it is more energetically expensive to 
have many of them so the local noise has to decrease. As 
for the exchange term, we hypothesise that it increases 
because the gravitational field creates a lensing effect 
such that more virtual particles emitted by detector 2 
hit detector 1. 

As we previously mentioned this effect scales linearly 
with the strength of the gravitational field M/ri, so for 
detectors with a ^ 1/AE and a > L p we have for 
the Earth Nm ^ 10 -9 -/V(o) while for the Sun we have 
-^(l) IO^A^o)- Since entanglement swapping from 
the entanglement of the vacuum has still not been 
observed, we expect that observing will be very dif- 
ficult. Nevertheless, it should be interesting to see if this 
effect can be modeled in a quantum field analog like a 
linear ion trap j!5j . 



VII. OUTLOOK 

Our calculations depended on the assumption that the 
vacuum of the quantum field is described by the Boul- 
ware vacuum. If we had considered two Unruh-DeWitt 
detectors near a black hole in an Unruh or a Kruskal vac- 
uum [8] , the Hawking temperature seen by both detectors 
would have increased the local noise significantly such 
that the entanglement between both detectors should be 
degraded, not enhanced. It should be interesting to inves- 
tigate in detail to what extent entanglement extraction 
by detectors near black holes and stars is affected by the 
properties of the corresponding vacuum states. 

In this context, it should also be interesting to inves- 
tigate whether one can effectively model a black hole by 



G 



using a confining potential, say on a shell. Indeed, a trap- 
ping potential can have horizons, so it may be possible to 
have a non-trivial vacuum in which particle production 
occurs because of the potential. Such an analysis could 
show the Hawking effect and its various open questions 
in a new light. 

Since we observed that the exchange term \X\ increases 
because of the gravitational field, it is tempting to spec- 
ulate on the Casimir-Polder force near a constant den- 
sity star. Indeed, the exchange term and the Casimir- 
Polder force have essentially the same interpretation, 
that is they are the result of a continuous exchange of 



virtual particles. We therefore conjecture that Casimir 
or Casimir-Polder forces can slightly increase in a weak 
gravitational field. 
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